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Abstract 



crystal) as well as the term introduced by Kardar, Parisi, and Zhang (KPZ) to 



■ The deposition dynamics of particles (or the growth of a rigid crystal) on 

ON 
OO . 

. a disordered substrate at a finite deposition rate is explored. We begin with 

o\ 
o 

■r^lj- ' an equation of motion which includes, in addition to the disorder, the periodic 

"j^ ■ potential due to the discrete size of the particles (or to the lattice structure of the 

a 

^ ■ account for the lateral growth at a finite growth rate. A generating functional 

> ■ 

"-H . for the correlation and response functions of this process is derived using the 

I approach of Martin, Sigga, and Rose. A consistent renormalized perturbation 

expansion to first order in the non-Gaussian couplings requires the calculation 
of diagrams up to three loops. To this order we show, for the first time for this 
class of models which violates the the fluctuation-dissipation theorem, that the 
theory is renormalizable. 
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We find that the effects of the periodic potential and the disorder decay 
on very large scales and asymptotically the KPZ term dominates the behavior. 
However, strong non-trivial crossover effects are found for large intermediate 
scales. 

I. INTRODUCTION 

1. General. 

A few years ago the near equilibrium dynamics of a growing crystalline surface was elu- 
cidated [0,0. It was found that a roughening transition occurs at T = between a high- 
temperature rough phase and a phase with a flat surface for T < Tj.. The mobility of the 
growing surface drops from a finite value to zero as T ^ T/. In the low-temperature phase the 
growth is "activated" with the formation of higher "islands" on top of the flat surface. Similar 
behavior occurs in deposition of cubic particles with diffusion on a fiat substrate. The transition 
in this case is as function of the inherent noise due to spatial and temporal fluctuations in the 
deposition. 

In view of the existence of a low temperature (or low noise) flat phase, the question of 
how disorder in the substrate will change the behavior had to be addressed. We have ini- 
tiated a comprehensive study of the related questions. A short letter which announces the 
novel and surprising results was published elsewhere 0. In a previous full paper [Q (denoted 
by I in the following) we have presented the detailed calculation and analysis for the near- 
equilibrium dynamics. In this regime the averaged growth rate is small. The equation of 
motion is derived from the Hamiltonian of the system. Detailed-balance and the Fluctuation- 
Dissipation Theorem (FDT) ||^ both hold. Yet we found very non-trivial results: super-rough 
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correlations, temperature-dependent dynamics exponent, and a non-linear relation between the 
average growth rate and the driving force, were found below a super-roughening transition 
temperature Tgr- 

In the present paper, the second in the series, we present our detailed calculations and 
results for the dynamics far-from equilibrium. In this case the equation-of-motion cannot be 
derived from a Hamiltonian. Detailed-balance and the FDT are both violated. That situation 
represents a much more serious theoretical challenge since even the very renormalizability of 
the process is questionable. 

The equation of motion we analyze describes the deposition of cubic particles on a random 
substrate. It will also apply to the surface of a crystal if the solid is very rigid. Since the effects 
of the disorder in the substrate will be felt only up to a height h* (which is larger if the solid 
is more rigid), our theory will apply as long as h < h*. 

The width of a growing surface w follows generally the scaling form [§,0: 

w{L,t)^L'^f{t/L^), (1) 

where t is the time and L is the linear size of the system, a is the roughening exponent and z 
is the dynamic exponent. 

The leading difference between near- and far-from equilibrium is due to the lateral growth. 
If the growth rate is finite the lateral growth add a term proportional to (V/i)^ to the equation 
of motion. This term was derived by Kardar, Parisi and Zhang (KPZ) as the most relevant 
term in the expansion in term of Vh. 

As the RG analysis shows, the KPZ nonlinearity is marginally relevant in 2+1 dimen- 
sions. The asymptotic behavior is controlled by this nonlinearity, which violates the FDT. 
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Consequently, the roughness of the growing surface is also affected by the nonlinearity. Several 
simulations showed a ~ 0.4 [9,10] which is larger than that of the near-equilibrium case (a = 0). 
Many variants of KPZ-related models have been studied in recent years [6,7,11-22]. 

The effects of such a term on the growing crystalline surface on a flat substrate was consid- 
ered by Hwa, Kardar, and Paczuski (HKP) |Q. We review and revise their results in the next 
section of this Introduction (Ch.I). Ch. II will be devoted to the derivation of the generating 
functional, and the description of the RG procedure. Ch. Ill, IV and V are devoted to the 
outline of the calculations of the different renormalization factors. Extensive details of these 
calculation are relegated to the appendices. In Ch. VI the recursion relations are derived. Ch. 
VII is devoted to the analysis of the asymptotic and crossover behaviors which follows from 
these recursion relations. Our main conclusions are summarized in the last chapter VIII. 

2. Review of previous works 

The system we shall study in this paper has three important and non-trivial ingredients: 

(i) The periodic potential, (ii) The disorder in the substrate, and (iii) The KPZ non-linearity 
arising from the lateral growth in the presence of a finite driving force (or deposition rate). 

Paper (I) discussed the case in which (i) + (ii) are present. The major physical consequences, 
namely the existence of a super-rough "glassy phase" for T < T^^ with intriguing static and 
dynamic properties, were discussed there and will not be elaborated further here. Our goal 
here is to see how the addition of (iii) (i.e. the KPZ term) modifies the behavior. 

However, we can also take another point of view and ask how the addition of (ii) (i.e. the 
disorder in the substrate) is modifying the behavior found in presence of (i) + (iii). The model 
which discussed the non-equilibrium growth on a fiat surface in presence of both the periodic 



potential and the KPZ non-linearity was analyzed by HKP 23 . The equation of motion which 



describes the growing surface under these calculation is: 




F + iy{V'h{x,t)) + ^{Vhr + ^ 



sm[^{h{x,t))] + C{x,t), 



(2) 



In this equation h{x, t) is the local height at time t, x — {x, y) are the coordinates in the 
2d basal plane; /i is the microscopic mobility; F is the driving force; u is the surface tension; A 
is the coefficient of the KPZ non-linearity (which is proportional to F); 7/2 is the coefficient of 
the leading harmonies (higher harmonies are irrelevant and the lattice spacing is taken to be 



spacing; and r]{x.t) is the noise in the deposition (or due to thermal fluctuations) which obeys: 



They obtained a set of recursion relations from which they reached several conclusions. 
They found a critical temperature Tc- For T > Tc, y2 decays slowly to zero and the large scale 
behavior is determined by the KPZ coupling. Approaching Tc from above the linear-response 
macroscopic mobility (Namely the ratio v/ F where v =< ^ > is the averaged growth rate in 
the hmit F ^ 0) vanishes as (In \T - Td)""- 

In the low temperature phase T <Tc they have found that 7/2 grows indefinitely large and 
therefore they concluded that the surface is flat. 

While reanalyzing their work we have discovered a term that was overlooked in their cal- 
culations and which might affect their latter conclusion. Indeed it may be shown that a term 
of the form yi cos(27r/i) is generated under renormalization from the contraction of the terms 
A(V/i)^ and t/2 sin(27r/i). This term also feeds back into the renormalization of ^2- 

The recursion relations to lowest order are: 



unity for simplicity); 7 



Yi where h is the vertical lattice spacing ;a is the horizontal lattice 



< r]{x,t)r]{x ,t) >= 2Dd{x- x)5{t-t) 



(3) 
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dyi ,^ Ti^iD X 

-rr = (2 )yi - p-2/2, (4) 

al V V 

dy2 ,^ Ti^iD X 

-7]-= 2 )y2 + p-yi, 5 

where p = |[47r^ ln(4/3)(^)^ + (^)]. The two harmonic terms may be combined into a single 
term: 

\y\sm[27rh + 'd{l)] (6) 

with = yf + and = tg~^yi/y2. 

If we look at the flow of \y\ and ?9 we find that indeed \y\ oo for T < Tc which would 
indicate a flat surface if i} was to remain constant. However, the recursion relations imply that 
the phase shift angle is rotating with / like: = wl with an angular velocity w ^ X/u. Since 
/ = In 6 where b is the rescaling factor that means an ever changing ■(9 ~ w In 6, with rescaling. 

Therefore it is not obvious that the surface is flat. Hopefully, higher-order terms in the 
recursion relations of yi,y2 and u will help to identify with more confldence the nature of 
low-temperature phase. 

3. Introducing the disorder: 

As we explained in I, and as is clear from Fig. 1, the disorder will shift the origin of the 
periodic potential by a random and uncorrelated amount at every point x in the basal plane 
ii- 

We take the deposited particle to have a rectangular shape with a square base of linear extent 
a and height b in the growth direction. Then the periodic potential becomes proportional to 
sin[^(/i(x, t) +d{x))], where d{x) are the local deviation in the height of the substrate. Let us 
denote the associate phase 0(x) = 27id{x)/b. It obeys: 

< e-®(^")e-®(^^) >= a^S\x-y). (7) 



We have assumed that d{x) are typically of order b (or larger) and that if correlations exist in 
d{x) at different x they are at most short range (in which case they fall in the same universality 
class as the 6 correlated disorder we study here) The equation of motion we need to study is 
therefore: 

^^,dh(^ = F + u{V'h{x, t)) + ^{Vhy + ^ sin[7/i(f , t) + e(x)] + C(f , t). (8) 

II. GENERATING FUNCTIONAL AND BASIC DIAGRAMS. 

The Martin, Sigga, and Rose (MSR) ||2^-p6| method is utilized to obtain the generating 



functional for the correlation and response theory. An auxiliary field ih{x, t) is introduced 
to enforce the equation of motion through an integral representation of the delta function. 
Then the thermal noise and the disorder are averaged upon to yield the following generating 
functional: 

{Ze[J, J])disorder = J VhVhexpiJ (f xdt[D fi^fi" - h{—h - fivV^h - li-{Vhf)] 
+^Y^ J J ^(^' ^)^(^' cos(7(/i(f, t) - h{x, t'))}. 

(9) 

a. The Basic diagrams: 

As explained in details in I the Gaussian (quadratic) part of the "action" gives rise to the 
bare response function 



u;)4>{-q, -Lj)) = (10) 



which is depicted in Fig. 2. 



The bare correlation function is : 



2' 



where m, the mass of the field 0, is introduced to control the infrared divergences in the Feynman 



integrals. (In the two-dimensional regularization, it is notorious |25| that their infrared and 
ultra-violet divergences will mingle together without introducing masses for the fields.) 
In the momentum and time representation, they are given by ||26|| : 



{(j){q,t)4){-q,t')) = e{t - t')e-M('?^+™')(*-*') 



(12) 



{mt)<P{-q.t')) 



q2 _|_ ^2 



-At(g2+m2)|t-t'| 



(13) 



where 9{t) = 1 for t > and 9{t) = for t < 0. The bare vertex (j)(j) cos(0 — (j)) is drawn in 
Fig. 3. 

b. The renormalization group procedure: 

We follow the minimal subtraction scheme. The renormalization parameters are related 
to the bare ones by the so-called Z factors. Minimal scheme consists in extracting from the 
diagrams only the divergent parts which are all expressed in terms of functions of the dimen- 
sionality of the system. 

The bare and renormalized vertex functions can be related by factors of Za,-, Zi. For instance, 



(14) 



where and (^o label renormalized parameters {g,^,---) and bare parameters {go, fiQ, ■ ■ 
respectively, q and u are the external momentum and frequency, respectively. In the corre- 
sponding vertex function, a is a short-distance cutoff, and k is a mass scale. Here T^^l stands 



for the vertex function with L external (f) hnes and external (j) lines. The factors, and Z^, 
are set to remove the divergent parts of the vertex function F. 

The following renormalization constants are defined through the relations between the bare 



and the renormalized couplings 1 25, 27 -29 



Dq = ZdD, go = Zgg, Xq = ZxX, (15) 

ml<p' = m'<Pl, 7oV = 7% = Z^<j)l, = {Z^)4>1, (16) 
where Z^ = Z^-. 

In the next chapter we concentrate on the procedures to calculate the Z factors. The details 
are relegated to the appendices. 

The renormalized perturbation theory even to lower non-trivial orders in g and A has to be 



consistent order by order in 7. It also requires the calculations of Feynman diagrams p5| , ^ 
with up to three loops. 



III. CALCULATIONS OF Zd AND 

The renormalization of n is not affected in the presence of the KPZ nonlinearity since the 
associated vertex function Ti i comes with one external (f) and one external (j) and the basic 
vertex A contains derivatives on its two (f) legs. Thus the factor remains the same as in 
paper I (equilibrium dynamics), and so does the recursion relation for /i. On the other hand, 
the parameter D suffers additional renormalization of order A^. Basically, the renormalization 
of D from A^ is the same as that of D encountered in the KPZ model. Here, we still focus on 
the same vertex function F2,o as we did in the previous paper (I). 
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Obviously, the first nontrivial contribution begins from tlie second order in A. As shown in 
Fig. 5, the vertex function is modified by the associated integral with A^. The corresponding 
integral is given by: 

AioXo ,J_ rco 1 ^2Dofil)Y-p^ 

1 /■°° 1 1 

= \xl{Dlf^l)[-l-ln{cm'a% (17) 

As illustrated in Fig. 5, the symmetry factor for this diagram is 2 and another factor 2 arises 
from the differentiation of the external legs. One more factor | is due to the coefficient from 
the expansion of interaction with power 2 for the utilization of A^. With the combination of 
Eq. (D6) in paper I, we have: 

- 2D^2 ^ (^,)2[_2Z^^^2 _ :^^2(_^i^^^^2^2)) ^^2y^^^^l^^^^2^2)]_ ^^g) 

Therefore 

DflX"^ 1 2 2nn , I'^Vcg . I 2 2n /in\ 

= 1 H — (— ln(cm a )) + — — — ln(cm a ). (19) 

4 47r 2iJ/i 



IV. THE CALCULATION OF Zg 



As in the harmonic model of paper (I), we consider the vertex function r2,o(^, ^; iii 



the limit \t — t \ oo |28]. The calculation of Zg will be based on Eq. (8.3)in GS |28], where the 
renormalization factor is substituted in that equation to make r2,o(if, ^; (<? ~^ 0) finite. 

Therefore to find out Zg is just to calculate the renormalization of r2,o- 

The contributions to Zg we need to sum are of order g"^ (as in paper I) and of order X^g. 
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The combination of g and leads to two types of diagrams, 2-loop and 3-loop diagrams. 
Some of the associated diagrams are canceled by each other as shown in Fig. 6 and Fig. 7. 
The other non-vanishing diagrams, including six 2-loop and two 3-loop diagrams, are shown in 
Fig. 8 - Fig. 15. The detailed calculations of 2-loop integrals are given in Appendix A, where 
we also explain the cancellation of sub-divergences of some diagrams. The sum of the leading 
and sub-leading divergences contributing to Zg are listed in Eq. |22| (see the third and fourth 
terms). 

In the Appendix B, we present the detailed calculation of the leading divergences in the 
3-loop diagrams and also show that they do not contain any sub-divergence. 

Now, what remains to complete the 3-loop results is just to sum up the leading divergent 
terms in In(cm^a^), which essentially contribute to the recursion relations. The In(cm^a^) 
contribution of the diagram in Fig. 14 will be 

-^[(3.1-1) -2(3.1 -2) -2(3.1 -3)] 

= hs + 2Io + 4/^ + 20l In(^) - Q{H^)f - 4 ln(2) - ^(In 2f + 22 In(^) 

+3$(l,2)-3<l>(i,2)-S(i,2) (20) 

The contribution from the diagram in Fig. 15 is: 



[-(3.2 - 1) + 2(3.2 - 2) + 2(3.2 - 3)] 

= ^[^^c + Id-2Ij, + \ In(^) - \ In(^) - \ ln(^)]. (21) 
Now we are in a position to calculate Zg. Inserting the above calculation results into the 



self energy in reference 28 , we obtain: 



y . ew 2 2^ , A^^(D/i)^ [In(cm^a^)]^ (-5 + llln(|)) 

Zg = l-5Q ln(cm a ) + ^7 {^lD) A 
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^"^'"^'"'^ + (-90.5)y(D/.)^A^M^. (22) 



(47r)2 ^ . / V (4^^ 



By using 60 = 6 + iA^(D/i)^7^^^^^^y|-^, we can calculate the recursion relation of g, as will be 
shown in the next section. As we will see the term ln(cm^ ay is canceled when we derive the 
Callan-Symanzik (CS) equation. Thus the scaling equation is consistent. 



V. THE CALCULATION OF Z 



A 



For the calculation of Zx, one should consider the vertex function Fi 2. For the diagrams 
in Fig. 16 and Fig. 17, one can write down the associated integrals as: 

nla = r £k r dn „ 2D/.V- (f -k) ^ 

and 

nib = f £k r dn . . 2D,^p.il-k)^ 

where 



00 
00 



— CXJ /i 

1 



(tQ) / dte'™(e^ - 1). (25) 



/i^D ' ' Jo 

For simplicity, let x = (| + /c)^, and y = (| — A;)^. The summation of nla in Eq. ^ and nib 

in Eq. ^ is proportional to '^^ (x2+n2Ki;^+n2) /(^)- 

With the help of Eq. (^), the frequency dependent part can be integrated out first: 

r dp. {x){i^)e'^^ _ e-^'* _ e-y^ ^^g^ 

J -00 (x^ + fi2^(y2 _j_ Q2-) rjp. _ y2 ^2 _ y2 
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We then have: 



^ / r)2 ,70,o\/— *7\o'- 



(27r)2y^oo 4(2j + A;2 + m2)(p. A?)/i2 
x/i[(f + ^r + m2] 

= p.{P- k){-2)- _ j^^e 

^4 ^ (p. A;) ^^2 ^ ^ + ^2 + ^2 ^ ^ 



(27) 



In the hydrodynamic (long wave-length) limit, p — > 0, The relevant term in the first term in 
Eq. ( p7[ ) can be easily found as: 

c;2^(_2)^-te-^('='+™')* = -p2— e~^"'*. (28) 

-00 2 2y(i 

The relevant term in second term of Eq. ( pTj) is: 

p2(2) / dOcos'e kdkt—- ^ =p2_g-/.m2t_ (29) 

JO Jo + 777.2 2/i 

Thus there are no contributions to the renormalization of A due to their mutual cancellation. 
Other possible diagrams arise, but result in no contributions. In Fig. 18. those two diagram 
will not contribute the renormalization when one impose the long-time prescription {Qext = 0). 
The diagrams in Fig. 19. do not contribute neither since the interaction g is local in space and 
therefore there is no j7-dependent part of the vertex ri^2- 

To sum up, A suffers no renormalization within the perturbative expansion to order of g, 
and therefore Zx = 1- 

VI. RECURSION RELATIONS 



Once the Z-factors are known to the leading order in g, the recursion relations are obtained 
via the so-called /3-functions [ ^p6| , p9| : 
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/3, = .(^). = -^?-(^). (32) 



Pu = ^{-^)b (33) 

where subscript b means that all bare parameters are fixed when one performs the differentia- 
tions [2^,26,2^ and k is a mass scale. The renormalization of the couplings may also be related 



to the same P functions. 

The renormalization Z factors are the ratios between the corresponding renormalized and 
bare parameters. Therefore it is a standard procedure to extract from their dependence on the 
momentum scale k, (or the bare mass mo) the fiow of the renormalized couplings under rescaling 
of all length scales by 6 = exp(/). The first step is to compute the so-called beta functions 
which when subtracted from the naive (engineering) dimension of the couplings yield the fiow 
equations. Following this procedure, rescaling length scales x — > bx, momenta k b~^k, time 
t — >■ b^t and frequencies w — >■ b~^w, we find the following recursion relations: 



du TT'-f 



2 



(35) 



dl 4z/(D/i)3- 
dF 

^ = 2F + nX (36) 
dl 

^ = i-D,+ '-^)D (37) 

-dl = ^''^ 
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7 is not renormalized because Z^ = 1 and keeps its bare value 7 = 27r. For the same reason 
V is not renormalized and may be chosen as v = 1. 

The two constants are c = \e^^ = 0.7931 where E is the Euler Constant, and c' ~ 180.08, 
which is derived from the sum of the terms contributing to Zg in Eq. (|22|) . 

VII. THE ASYMPTOTIC BEHAVIOR 

In this chapter we proceed with the analysis of the physical implications of the recursion 
relations. We begin, in the next section by looking at the asymptotic destination of the flows 
which will yield the physical properties on very large scales of time and space. In the following 
subsection will analyze the crossover behavior which determines the properties on intermediate 
scales. 

1. The asymptotic behavior 

The analysis of the recursion relations may be facilitated by the introduction of a "tem- 
perature" like variable (temperature is not well-defined far away from equilibrium where the 
Einstein relation does not hold). Here we define it by T = Dfi (it is not the thermodynamic 
temperature). Its recursion relation is obtained from Eq. |3^ and Eq. |38|.It obeys: 

where the critical value Dfi = l/vr is substituted. Since A and 7 are kept constant, this equation 
may be integrated: 
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T{1) = Toe^'. (42) 
We see that no matter how small Tq is T{1) will grow indefinitely with I — Inb such that: 

T{l)=To{-f/'^\ (43) 

So the efi^ective "temperature" becomes higher on longer length scales. Asymptotically the 
system is always at a high "temperature". The growth of T is, however, quite slow. Therefore 
crossover effects discussed below plays an important role. 

What is the effect of high T? For that we have to look at the flow of g: 

It is clear that if T{1) grows very large it will cause g{l) to decay to zero, no matter what 
are the bare values go, Tq and A. Once g ^ the asymptotic behavior becomes equivalent to 
that of the KPZ equation. 

We, thus, conclude that asymptotically on very large scales and very long time the scaling 
properties are these associated with a far-from equilibrium growth without the disorder and the 
periodic potential. The behavior will be determined by the effect of the lateral growth alone. 
The KPZ properties in 2 + 1 dimensions (dominated by an inaccessible fixed point) will be the 
asymptotic ones for the system under consideration. 

The Crossover Behavior As we have found in the previous section, the "temperature" T{1) 
rises with the scale quite slowly. Hence the decay of g to zero might also be slow. As long as 
g is not vanishing the effects of the disorder and the periodic potential are still felt. Hence we 
should expect a slowing down of the dynamics. This slowing down will be observable on larger 
scales as well because the mobility obeys the equation: 
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dl T{1) 
and therefore 



(45) 



^(/) = /ioe-^^v^/o'^''[^('')/^('')l (46) 
. The ratio yu(/)//io does depend on the integrah 

m-jl^/. (47) 

and clearly J {I) is sensitive to g{l) on small scales as well. 

To evaluate J{1) we need to know T(l) given in Eq. ^ and g{T) which we calculate next. 
Given r(/), g{r) is found by integration of its recursion relation: 



^7(0 ^7(0) 87r2 





where 

s{x) = - 2]a; + To^(e^'^'/2^' - 1). (49) 
Clearly the second term dominates for large x. Hence we see that for large /: 

c/(/) ~ exp{-exp(A/)} (50) 

with A ~ ^. 

It is easy to see that for large enough I, g{l) decays to zero faster than exponential. Since T{1) 
diverges, the most important contribution to J{1) comes from small (or at most intermediate) 
values of /. At large /, J{1) approaches asymptotically a constant and its dependence on / 
becomes much weaker. 

This asymptotic value of J{1) depends mostly on the bare values of the parameters. To 
summarize, the mobility decays fast on initial small scales and then saturates to an almost 
constant value on large scales. 
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VIII. CONCLUSIONS 



In this work we have investigated the behavior of a class of growth systems in which three dif- 
ferent effects play important roles: Periodicity, disorder and lateral growth. Our model applies 
to the situations in which all three effects are present in deposition processes or solidification 
of rigid crystals. 

In I we looked at the system near-equilibrium when the lateral growth is negligible. There 
we found a continuous transition from a rough phase at high temperature into a super-rough 
and glassy phase for T < Tg. 

The main conclusion of the present work is that far-from equilibrium the KPZ term prevent 
this transition. The ultimate asymptotic behavior is dominated by the KPZ non-linearity while 
the second non-linear term, (obtained upon averaging the periodic potential over the disorder) 
is irrelevant. 

We have seen, however, that while this term is decaying it still affects the behavior on 
intermediate scales. The effect of the disorder in the periodic potential is to slow the dynamics. 
In particular the mobility decays from its bare value as: 

^ = exp[-47r2xl.78J(/)], (51) 

where J(Z) is given by Eq. 

It is clear that the behavior on intermediate scales drastically depends on the bare values 
of the parameters. 

We have also shown that the theory is consistently renormalizable to first order in A and g. 
To obtain the correct renormalization we had to keep diagrams up to and including three non- 
trivial loops. As far as we are aware of no other calculation have shown the renormalizability 



up to this order for a dynamic system for which the fluctuation-dissipation theorem is not 
satisfied. It is reassuring to see that the renormahzation group can be successfully applied to 
the dynamics far-from equilibrium. 
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APPENDIX A: 2-LOOP CALCULATION FOR Zg 



As explained in Chapter IV, the calculation of Zg is based on the evaluation of the vertex 
function r2,o- 

For clarity, we neglect the prefactors and symmetry factors. Here we have 6 2-loop and 
2 3-loop Feynman diagrams. The 2-loop diagrams will be denoted by FD 21#, and 3-loop 
diagrams will be denoted by FD 31^^, where ^ stands for the sequel. First we look at the 2- 
loop diagrams. The basic rules for the calculations of these diagrams have been described 



in the Appendix of our previous paper I. Here we shall simply write down the corresponding 
integral for each diagram. 

We employ the momentum-time representation for correlation and response functions, in 
terms of which the corresponding integrals will be easily handled. 

Feynman diagram (FD) 211 is shown in Fig. 8. The corresponding integral over time is 
given by: 



FD 211 



^Jo ^ {q^ + rn?)^^^ + m?) 



The integration of the time dependent sectors gives: 



^io 



(Al) 



y 



^ [g^ + (p — g)^ + — p2] 



2(p^ + m?)[p'^ + g^ + {p — q)"^ + 3m^ 
By decomposing p- q into 
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(A2) 



p- q = — + + {p- q)^ + Sm^ - 2p'^ - 2q^ - Sm^] (A3) 
2 



and q- (q — p) into 



(g^ - p- q) = ^[p^ + q'^ + {p- g)^ + 3m^ - 2p^ - 3m^] (A4) 



we obtain 

{q-p) 



{p^ + m?){q^ + m2)[p2 + g^ + (p — g)^ + 3m2] 

1.1 2 



2 (p^ + m?){q'^ + m^) (p^ + m?){q^ + m2)[p^ + g^ + (p — g)^ + ?>m?\ 
2 

irrelevant terms]. (A5) 



(p^ + m^)^[p^ + + (p — g)^ + 3m^] 

Since we impose the minimal subtraction (MS) P3| , |2U| scheme, the unwanted nonsingular 
parts (finite parts) will be ignored in all calculations in this paper. In this whole paper, an 
irrelevant term means a term which does not contribute to the singular part of the integral. 
We sometimes use equal sign to represent the equality of the singular parts on both sides. 

We substitute Eq. (|A|) and Eq. (|X5D into Eq. ([Al]), and obtain: 

Eq. (Al) 

_ 1 ^ {(f - P - Q) [p^ + q^ + {p- qf + 3m^ - 2p^ - 3m^] 



4 {p^ + m?Y{q^ + m?) {p^ + m?){q^ + m?)[p^ + g^ + {p — g)^ + 3m^] 
[p^ + g^ + (p — g)^ + 3m^ — 2p^ — 3m^^ 



(p2 -f m2)2[p2 + g2 _|_ g*)2 _j_ 377^2] 



-} 



Ir 1 P-Q 



+ 



4 (p^ + m^)^ (p2 + m2)(g^ + m^) (p^ + m^)(g^ + m^) 
2 1 



(g^ + m2)[p2 + g^ + (p — g)^ + 3m^] (p^ + m^)^ 
2 



(p2 + m^) [p2 + g2 + (p — g)2 -)_ 3777,^] 
4^ 



^-B-Cl (A6) 



where 



21 



(p2 + m^) [p^ + + (p — g)^ + 3m^] ' ^ 
and the second term in Eq. (|A6|) is discarded because it vanishes after the integration over the 
momentum variables. 

In the same way, the associated integral for Fig. 9 reads: 

FD 212 = rdt, dt [g-(p-g)] ^-(WH. 

To begin with, we integrate over the time variables, and that yields: 

J-oo ^ 



JO 



{p^ + rn?){p'^ + + {p — qf' + 3m^) 
By inserting Eq. ([A10| ) into Eq. (^), we obtain: 



[q^ + m?){p'^ + m^) [(p — g)^ + m^] [p^ + + (p — g]^ + Sm^) 

_ 3 

4 + m2)[(p — g)2 + m2][p2 + g^ + (p — g)^ + Sm^] 
1 1 



2 [p^ + g^ + (p — g)^ + 3m2](g2 + m^) 



where 



(AlO) 



(All) 



^ = 1 (A12) 

(p2 + m?)[{p — qf' + m2][p2 + g2 _|_ _|_ 3777,2] 

For the Feynman diagram in Fig. 10, the integral with the time variables is written as: 
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iv^ + rri^^ia^ + m?) Jo Jo 



(p2 + rn?) (g2 + m?) Jo Jo 
Again, the time dependent sectors are integrated out in advance. 



Jo ^ 



2 (p^ + m^)[p^ + q"^ + {p — q)"^ + Sm"^] 
With the help of Eq. (p^il) , Eq. (|XT3D is simplified into: 



(A14) 



{p-q){p-q) 1 1 



(p^ + m^) (g^ + m^) 2 (p^ + m?) [p'^ + + (p — g)^ + "im?] 

= ][^C-B]. (A15) 

Now we turn to the diagram in Fig. 11. The associated integral is represented by: 

FD 214 = HdU r rft,e-(^^--^)'^ 



(g^ + m^)[(p — g)^ + m^] jo ^ 

g-[gH(p-g)2+2m2]|t^|g-(p2+m2)fe-tj;)_ (A16) 

The time dependent sectors in Fig. 11 is identical to Eq. (|A10|) , so we will not repeat the 
calculation here. 

In the same fashion, the integrand takes the form: 

(g2 + m^) [{p — q)"^ + m?] [p^ + g^ + (p — g)^ + ?>rn?] 

= [\b-A]. (A17) 
For the diagram in Fig. 12, we have the associated integral as: 
FD 215 



{q^ + m?)(jp' + m?)Jo ^Jo ^ 

(A18) 
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The time dependent term is integrated out first: 



Jo Jo ^ 



Then we obtain: 



(A19) 



1 1 -\b-2C). (A20) 



(g2 + rn?){p'^ + m^) 2 (p^ + rn?)\p'^ + + (p — g)^ + Sm^] 4 

We turn to the simplest figure among 2-loop diagrams, Fig. 13, of which the related integral 
is evaluated as: 

/ dtp'Ro{p,t)Co{p,t)y = -[ d'p (A21) 
-oo Jo 4 J-oo (p^ + m^) 



1 roo 1 POO 

FD 216 = - - / d'^p— / d'^p- 



4 7-0O (p^ + m^) J-oo (p^ + m^) 
= -^S (A22) 

Now we go on to perform the integration over the momentum variables. Before the evalu- 
ation of the singular parts of integrals, it will be helpful to present some identities, which will 
play important roles in later calculations and have been frequently employed in these types of 
calculations. The first one is the Feynman parameterization formula, which reads: 
1 r{a + p + -f + --- + e) n ri 



X 



{Ax + By+--- + £;z)"+/3+7+-+^ • ^ > 

A set of integral formulas is also valuable and is shown as below: 

Jo = / d'^k^ ^ = 7^(M - p'Y'^-'^Via - ^). (A24) 
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d'^k^ ^ = -pVo. (A25) 

(A;2 + 2A;-p + M)« 



oo 



With these formulas, one can evaluate the singular parts of the integrals. Let X — J J d'^pd'^qX, 
where X = A, B or C, and d = 2 + e. The evaluations of B, C and A are carried out as below: 



/■CO 1 



oo (p2 + m^) 

= 7r^/^(m^)f-^^^ (A26) 

/oo ^ ;>oo ^ 
-oo ^ J-oo ^(p2 + m2)[p2 + g2 + (p — g)2 + Sm^] 

= n ^'^P ^'^^Tl~^ ^ / 2 I 2 ^ ^ I 3 2^ + non-singular terms 

^ 1 r 1 rf/2 r(l-ci/2) 1 

2 7-00 ^(p2 + m2)^ r(l) (fp2 + 

_ 1 , r(-e/2)r(-6) 3 ,2^2e .A27) 
"2^' r(l-e/2) ^4^ ^^^^^ 



oo /-oo 



A= dy d^q 



oo i-oo -I- m2)[(p — ^2 _|_ 777,2] |p2 -)_ gr2 _|_ (p*— g^2 _|_ 3j^2j 

/oo /"OO 
rfV / d'^q- 
-oo J —oo 



ijp' + m2)[(p — g)2 + m'^]\p'^ + q'^ + {p — q)"^ + 2m2] 

1 /-oo /-oo /•! pi 

- d'^p d'^q / dxdy 

2 J-oo ^-oo ^0 Jo x+y<l 

mi 

{(1 — X — y){p^ + rn?) + x[{p — qf' + m^] + y{^p^ + — p - q + m^)}^ 

-oo ^ ^i-oo ^ Vo ./o x+y<l^^^^ [P^ - P • ?t2x + + (x + y)g'2 + m^]^ 



oAwi r(3) {Hx + y)-(x+i)V + m^''-'-'^ 



r(3) Jo -'0 x+v<\ J-oo 



r(3) Jo Jo x+y<l " J-oo " [{X + y)-{x+ |)2]2-^/2 
1 



X 



/o2 J inz 1.2-e/2 

,r(3-rf/2) /•! /•! 



,/2U3-d^ r i 

^""^ r(3) io io ^ [{x + y)-{x + f )2]2-/2 ^ 

I 1 d-t-2-2(2-e/2) / ^(1 + d/2) T{-e) 

^ [(x + i/)-(x + f)2]V2-r ^ ^ r(d/2) r(2-e/2) 

(^)d^^!l±ll31^iz5lf' [' dxdym'^- ^ UT^^T^- (A28) 

^ ^ r(l + e/2) 2 JoJox+y<i ^ [{x + y)-{x+lff+'l^ ^ ' 



25 



The singular parts of integral in Eq. ( |A28|) can be evaluated by the changes of variables, 
X = st, y = s(l — t). Let A' denote the x, y dependent part in Eq. (|A28|) . We obtain 



"1 rl 

A' = I / dxdy 



Jo .+y<i [{x + y) - {x + m^+^/^ 



1 rl s 



dsdt 



Jo [S~ - |)]2+^/2 

1 fl I 

dsdt- 



Jo sl+^/2[l - (1 - |)s]2+^/2 

= lo lo '^'^;i^^ [l-(l-H]2W2 - 1} + r ds^^. (A29) 

/^From the form of Eq. ( |A29| ), it is not hard to see that the singular term as e — >■ is the 
manifestation of the singular behavior of the pole s = in the integral. To simplify the 
expressions of equations, the term 1 — | is symbolized by a. Eq. (|A29|) can be rewritten as: 

Eq. (A29) = £ I' ,.*_J^if!£Z_f^ + I (A30) 

The analysis of two terms in Eq. ( |A30| ) will be carried out in order. The first term in 
Eq. ( [A3 01 ) is recast into: 



1 /-^ 1 a /-i 1 a 

dsdt — — , „ , ,„ + / / dsdt- 



Jo (1 - as)2W2 Jq Jo (1 - asY+^Z' 

1 '•I , a f\ , Ms)]a 

dsdt- —777 + e / / dsdt- ^ ^ 



n 



JO 



l-as)2+^/2 Jo Jo (l-as)2W2 



1 rl 



a 



The second term in Eq. ( |A31|) is of order e and therefore is discarded. The first term in 
Eq. ([A3 ID is evaluated as: 



1 a 1 1 

1 /•! 1 

^{ / dt r2 r - 1}, (A32) 

l + e/2Vo (t-^)iW2 J' ^ ^ 



where 
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/■II /-111 1 

Jo '^\t-^^Y+^l2 = Jo "^^^[(1-1)1+^/2 " ^] + io "^^t^ 



1 | + (l-|)fln(l-|) + --- , ^1,^ 1 



111 1 

^^4tl+e/2 (1 _ |)l+e/2 + ^(^) + Z^/^ 

l„(i) - ! + 0(e). (A33) 



The third term in Eq. ( |A31|) is: 



1 '•I a , a 



f rfi Ml -(1-^)1+ 0(e) 



= 31n(^) + 0(e) (A34) 
The second term in Eq. (|A30|) equals to — ^. 

With the combination of the prefactors in Eq. ( [A28| ) and the results obtained above, the 
singular part of A is given by: 



A = (nn— - - In(-) - + -]+ finite terms , (A35) 

e 3 e e 

where 7 is the Euler number. The finite part of the subdivergence diagram can be neglected, 



if the scale equation is well defined. One can always scale it away. ||2^ In the short distance 
cutoff p5| , |28[] , they appears as: 



A = -i_[l(ln(cm2a2))2 - In(^) \n{cm^a^) + ^ \n{cm^a^)] (A36) 
i? = ^[(ln(cmV))2] (A37) 
C = ^/-{Hcm\')) + \ In(^) In(cm^a^)] (A38) 

Now we retrieve the symmetry factors for each 2- loop diagrams (see Table 
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TABLES 



TABLE I. Symmetry factors and integrals of 2-loop diagrams 



Feynman diagrams 


Symmetry factors 


Integrals 


FD.211 


16 


(iB-C) 


FD.212 


8 


(fA-iC) 


FD.213 


32 


(C-iB) 


FD.214 


16 


(iB-A) 


FD.215 


32 


(iB-^C) 


FD.216 


16 


-iB 
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As one can see from Fig. 8 - Fig. 13, only Fig. 9, Fig. 12, and Fig. 13 contain sub-divergent 
diagrams. One also can verify this from the results listed in Table 0. The leading terms in the 
diagrams without sub-divergences, such as Fig. 8, Fig. 10, and Fig. 11, are of order In(cm^a^). 
On the other hand, the leading terms of the diagrams with subdivergence as mentioned above 
are of order [In(cm^a^)]^. Furthermore, the leading divergences of the diagrams in Fig. 12 and 
Fig. 13 are canceled out by each other. They have the same type of sub-divergences, (see the 
subdiagrams enclosed by the boxes in their own figures) which are not present in the lower 
order (1-loop) calculation. Another diagram in Fig. 9 includes the sub-divergent diagram (see 
the subdiagram enclosed by the box in Fig. 9) which occurs in the 1-loop calculation for D. 
As usual, it will be canceled when one calculates the recursion relations, eventhough the Zg 
factors contain some terms like [ln(cm^)a^]^. As we will show in the next Appendix, the 3-loop 
diagrams do not consist of any sub-divergent diagrams. Thus, the subdivergence in the present 
expansion only occurs in 2-loop diagrams. The cancellation of sub- divergences in Fig. 12 and 
Fig. 13, and that of Fig. 9 and Fig. 3 ensure the renormalizability of this theory (at least 
at this order). The sum of leading divergences and sub-leading divergences are summed up to 
contribute to Zg (see the third and fourth terms in Eq. 

APPENDIX B: 3-LOOP CALCULATION FOR Zg 

In this Appendix, we shall present the calculations for 3-loop diagrams, which are mentioned 
in Chapter IV. 

The diagram shown in Fig. 14 representing the integral is given by: 
FD 311 
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-oo J ~oo J —oo I 



-[p- A;) ■ k] 



-oo ^-oo J-oo (^q2 jY^2^\^(j^ ^ + m'^]{k'^ + rn?) 

coo ft-u 



Jo ^ Jo 

(Bl) 

To simplify the calculation, we denote a — {q^ + m^), 6 = [ip ~ + ^^]) c — {k!^ + m^), 
d = [{p — ky + m^], and e = (p^ + m^). Along the same hne as our preceding calculations of 
2-loop diagrams, we integrate out the time dependent term first. 



Jo ^ Jo ' 

= / c^t^e-("+''+'=+'^)*ne-("-"-'')*^ - 11 

Jo " 'b + c-e 

fOO 1 

Jo ^b + c-a 
1 



(c + d + e)(a + 6 + c + (i)' 

Therefore 

Eq. (Bl) 

jd^ r jdt -[p-q\[p-k]-k 



-OO J —oo J —oo 



-oo 

1 



(c + d + e)(a + 6 + c + d) 

-1 /"OO r'X) roo . in . Ft 



I roo roo roo 

(— ) / d'^p / / d'^k 

2i J —oo J —oo J —oo 



2 J-oo J-oo J-oo (^gr2 _|_ j^2^^^2 _|_ j^2^|'^p*_ ^-^2 _|_ j^; 
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1 



+ ^2) + (^2 + ^2) + [(^_ ^2 + ^2] + ^)2 + ^2]} 

2p- q 

{q^ + m^)[(p — ky + m^]{(A;^ + m^) + (p^ + m^) + [(p — A;)^ + m^]} 

1 

{{q^ + m^) + (A;^ + w?) + [(p — + m?] + [(p — A;)^ + w?]} 

2p ■ 



{q^ + rn?){k'^ + m^)[(p^ + m^) + (A;^ + rn?) + (p — A;)^ + rn?] 
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X 



(B2) 



, (B3) 

{(g^ + m?) + (k"^ + m?) + [{p — q)"^ + m"^] + [{p — k)"^ + m?]} 

where we have used the identities below to simphfy the expression of the equation: 



c+d+e 2 c+d+e a+b+c+d 

where s = p — k. To simphfy the calculation, we treat Eq. (B3) as a linear combination of 3 
integrals, 3.1-1, 3.1-2, and 3.1-3. Namely, 

Eq. (53) = —[(3 - 1.1) - 2(3 - 1.2) - 2(3 - 1.3)]. (B5) 

The term denoted by 3.1 — 1 yields: 

3.1 - 1 



d'p d'q d'k- P ^ 

-OO J — oo J — oo I 



-OO J-oo (^g2 _|_ j^2^^^2 _|_ ^2^|^^p'_ ^^2 _|_ j^2j 

1 

{(g2 + m2) + (P + m2) + [{p- q)^ + m^] + [{p- kf + m?]} 

OO J-oo 2 i(lji(lji(lj JO Jo x+y<l J-oo 
Pj_Q_ 

[k^ + k ■ p{—2x — y) + p'^{x + y) + p- q{—y) + yq^ + m?]'^ 

1 {6} 2 J-oo J-oo JO JO x+y<l 
PjJ_ 

{[(x + y) — {x + |)^]p^ + p ■ <l{—y) + yq^ + m?Y~''/'^ 

(^)d/2 r(2-6/2) .1 .1 ^^^^ 

2 J-oo J-oo JO JO x+y<l 

p-q 1 

(g2 + m2)?/2-e/2 [^2 _ ^. ^+ |p2 + i^2]2-e/2 

(^)d/2 r(2-6/2) .1 ri ri ^^^^^^ 

2 J-oo J-oo JO JO JO x+y<l 

r(3-e/2) 2^-^/2 ^.^ 

r(l)r(2 - e/2) ?/2-e/2[^2 _ ^+ z^p2 + £^2]3-e/2 

r(2 - e/2)r(3 - 6/2) ^^-/2 f ■ p2 r(2 

V / -n/o _ /r,\ ,.5-^/9 rrzA /^\9i„9 , z„.9T9_^v / i 



r(2 - e/2) y2-e/2 _ (|)2]p2 + ^^2|2-e V ^ r(3 - e/2) 
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2X2 Jo Jo Jo "^^"^ i-oo ^y2-./2(n)2-e [^2 + H^]2-. 

^- / / dxdydz— — — — — m( — )^/^f'(7r)'*/^ 

4 Jo Jo Jo x+y<l y2-e/2n2-.L V^^^ J V ) 

r(2 + e/2)r(-3/2e) 
r(l + e/2) r(2-e) 

^r(2 + ./2)r( 3/2.) f ; 



(B6) 



where A = (a; + y) — (x — and □ = |A — 

The working principles of extracting the singularity of the integral is based on the separation 
of the singular contributions from different points. All the calculation here follow this scenario. 
However, one should be able to keep track of those highly nested procedures. In the last line of 
Eq. ( |B6| ) , the prefactors before the integral contains a leading singular-pole of order i, and thus 
we should extract the contribution up to the zero order in e from this integral. The extraction 
of the poles and finite parts of this integral proceeds as below: 

1 /•! /•! ^2+6 

dxdydz- 



Jo Jo x+y<l ^ |/2+.n2+./2 

1 /•! /•! r<:/2 



'^'Jo Jo .^yj'''^\^n[{x + y)-{x + If] - f 
Rewriting x,y as y = st,x = s{l — t), we inherit the simplified equation: 

Eq. (37) = I dz^ I ^^^^ ,e/2,e/2[,_,2(l_|)2_|,,]2W2 

' t t ^^^^i^ -,(1- 1)2 _|]2W2 - 



JO JO 



Recast Eq. (^) into: 



rl rl rl ^^72 I 

Eq. (58) = dz / dsdt—- ,^ ttt;^, (B9) 

^ ^ Jo Jo Jo W2 (1 - ^)2+^/2sl+f(l -5^)2+^/2' ^ ^ 
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where a = ^rr^^- Consider the integration over s first, 

1 



Jo 



ds 



1 sa\2- sa] + [1 - sd^f ln(l - sa) H fi 1 

r/s — 1- / fh — 

_ +io 



A^ + A2. (BIO) 



Ai can be represented as the sum of An and A12. 



Ai = An + A12 



Furthermore, the term An can be decomposed into the following: 



^11 = ^111 + ^112 

A112 should not concern us because it contains terms of at least first order in e. We only 
calculate Am as: 



A 



111 



-(l + e/2)(l-sa)iW2Jo 



Substituting Am into Eq. (POf), we obtain: 



.MwJ - (B13) 



1 /-l ^e/2 ]^ 

"^^JO + e/2)t^/2(l - |)2+e/2 [(1 _ aY+e/2 " 



/' /' 



'(l + e/2)W2(l-f)2+e/2[l_il^]l+e/2 



1 /-l ^e/2 



— / (iz / (it 



/o 7o W2(l+e/2)(l-f)2+^/2 
-Sun — -B1112, (B14) 
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where we separate the integrand into two parts: 



1 rl 



~ l + e/2 Jo Jo - f - I - 1)1+^/2 

^ ^'dz z^/'Cdt ^^ , (B15) 



l + e/2 Jo Jo (1 - - f - 1)1+^/2 

The integral over the variable t in Eq. ( pi5|) is given by: 



,1 ]^ ^1 



t^{l - f - |)1W2(1 _ I) 7o - f - 1)^+^/2 

= Cnii + Cni2. (B16) 

The first term Cim is substituted into Eq. ( [B15| ), and the contribution is denoted by -Bum, 
which reads: 

Bum = I I dtdz- ; — \—, -t- + ■ ■ ■ irrelevant terms 

^0 Jo (1-f -i)(l-T) 

= 4[91n(^) +41n2 - ^(ln2)^ + $(1,2) - $(1,2) - E{\,2)] (B17) 

The evaluation of the integral is quite straightforward although tedious. One can find the basic 
integrals in Appendix C, whose compositions will be used to represent those complex integrals 
encountered in Eq. (pi5|) . 



Here we take the evaluation of -Bum as an example and set aside the rest of similar calcu- 
lations. 

dzdt -z\- -rr — 5—1 75— 

- - 4 (l-f)J(l-t+f) 

1 /•! 1 (i_i) 1 

dzdt — ^ ^ ^ 



Jo A^(l-t-l) (1-f)^ 
^Jo U 4^ (1-1)2^ 4^ (l-iW 



(B18) 
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where A = 1 - t - ^. Let u = I - \. Eq. ( |B18| ) turns into: 



Eq. (518) = -2 X l/2[ / du ""^^ , + 



1/2 ^ ln(i + f ) , ln(i + f ) ln(| + |) 



ln(l + f) ln(l + f) ln(l + |) ^ ln(i + f ) 

M U 1 — U 



Each term in Eq. ( pi8|) can be easily represented in terms of the basic integrals listed in 



Appendix C. For example, 



4 

■ 1 I ^^^ 



21n(2) + 31n(-) (B20) 
1-1/2 Infi + 1 1 

/ du^^-^ = [2 ln2(2) - - ln2(2) + $(-, 2) - $(1, 2)] (B21) 
Ji u 11 



/■1/2 infi + 1 
/ du ^^^^ =ln^(2) + 'l'(-,2)-$(l,2) (B23) 
Ji u 2 



1/2 Infi + Q 

rfM = 2 ln(2) - 3 In(-) (B24) 



1/2 infi + 1 

<i„(l±M = _1 1„2(2) - 4.(1. 2) + 4.(1 2) (B26) 

Ji M 2 2 

1/2 fl + ■?/■) 1 

du^-—^ = ^-,l)-^l,2), (B27) 
M 2 

where $ and S are defined in the Appendix C. 

Now we turn to the calculation of -B11112, which is defined as -B11112 = Jo dtCnu- 

1 rl 



-B11112 = / / dzdt 



JO 



where the t-dependent integral can be separated into: 

1 , 1 
dt = 
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1 [l-(l-/3t)^+-/^] 

t^+' (1 - ^)l+^/2 
1 , 1 

dt- 



dt—— 





/5 



(1 - /5)i+^/2 
D1 + D2 + 0(e) 



(B29) 



with /? being The term Di in Eq. ( |B29| ) is 



dt- 



lo (1 - pt) 
After substituting the above equation into Eq. 



ln(l-/3). 



one has: 



1 z^/Hn{l-^J fi [ln(4 - z) - ln(3 - ^)] 



"^^^ (1 - 1)1+^/2 



dz- 



(1 



.4. 



J^-JB + 41n^(-). 



0(e) 



(B30) 



(B31) 



Inserting 1^2 in Eq. ( P29|) into Eq. ( P28|) , we have: 

1 , 2^/2 



(it- , 

t^+' 



1 



-1 /•! 

e ^Jo (1 - t 
1 4 



6 , In(l-t) 



- / dz- , 

2 A (1 - j: 



^0 (1 - 4) 



-j)|41n(j)-|/. + ^/d 



(B32) 



Now we move on to calculate the contribution of A2 in Eq. ( [B1U|) . By inserting into 
Eq. (PSP, one has: 



= [ [ [ dzdsdt—^—;- + f [ [ dzdtds 

Jo Jo Jo 



'gl+e^l _ ' 7q 7o 7o 



4 i i e 

41n(-)— + / / dzdti—)i- 
6 — e JO Jo — e / 

-ln(l-f) + ln(z)-ln(t) 
[ (l-f)2 



(B33) 
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Here we decompose the second term in Eq. (|B33| ) into three parts, Pi , P2, and P3. 



"1 ri ln{l - ^ 
/o Jo (1 - f 



Pi = / / d^dt- 



4 



,.,(-4 /^(1 -f) |i I /s,., (-1) 



-[4/^ + /A]-41n(^) (B34) 



Jo (1-1)2 
1 ln(z) 



(1-t) 
1 '•I ln(t) 
Jo (1^) 



Ic (B35) 



dtdz ^^ = (B36) 



Now we go back to finish the calculation of A12 in Eq. (|B11|) . 



.1 .1 ln(l - 
^12 = / / dzdt- 



Jo '^^^ (1 - 



'0 io ''^'''(l-f)^ X (1-1)2 

4^^ , . , ^S.4ln(l-f-|) 



4 /"^ 4 

[4/^ + /A + 41n(-)] + Jc+ / dt- 
6 Jo t 



t (1-f 



zt\ 10 



lo Jo t(l-f)(l-f 
[4/z5 + /a + 4 ln(^)] + Jc + Li - Ml + M2, (B37) 



where Li, Mi and M2 are analyzed below. 



= £ dtj 1,1(5 _ * ) + 1-4 ln(?) + - 47c|. (B38) 
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The first term in Eq. (|B38| ) is divergent. As we will show later, it will be canceled by a term in 



r,„|Mi±i) + !!^l_8/"*li„(i 

JX u Jo t 



= 8[2 ln(2) - 3 In(^) + ln2(2) + <l>(i 2) - $(1, 2)] - 8/^, (B39) 

where A = (1 — and u = 1 — |. 

Ma = / dt-\n( 1)(— + 1) - / dt - n( ) 

Jo t ^4 4 A ' Jo t M 4^ 

= 4/ ^-(^)M4 + 2)-yo^^t^^^4-4) 

= 4[2 ln(2) + 3 In(^) + 2 ln2(2) - \ + $(1 2) - $(1, 2)] 

-/;.tlln(^4), (B40) 

where, as mentioned above, the first term in Eq. (P38| ) is canceled by the last term in Eq. (|B40| ) . 
The term -B1112 in Eq. ( |B14| ) is calculated below. 



,1 „l 

Bni2 = — I I dzdt- 



l + e/2 7o Jo t^/2(l-f)2W2 
1 /-i . 1 



lo Jo (l-f)2 
41n(^)+0(e). (B41) 



The second term 3.1 — 2 in Eq. (|B3|) is calculated as follows: 



3.1 - 2 

/oo roo roo _^ 
d'^p d'^q d'^k 
-00 J-00 J-00 (g,2_|_^2^J|p'_^'j2_j_^2j|^^2_|_^2^_j_|'^2_|_^_^^2_j_2^2jj 

1 



p-q 



{[g2 + (p - g)2 + 2m2] + [F + (p - A;)2 + 2m2]} 
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r(i + 1 + 1) r r r ?JL^ /' /■ ,,,,, 



4[r(l)]3 J-oo y~oo J-oo {q'^ + m?)Ji) Jo x+y<l 

{{1 — X — y) [(p — /c)^ + m^] + + — k ■ p + m^] 



+y{p +q +k — p-q—p-k + m)} + irrelevant terms 

^ r d'p r d'qjij^, f f dxdy f dzi^r/^ 

4 J-oo J-oo (g^ + m^jJO Jo x+y<l JO 

r(3 - 1 - e/2) 1 

f(3) {[1 - (1 - + ^. ^(_^) + yq2 + ^2}2-e/2 

^ ^ 4 i-oo ^io Jo x+y<i (g2 + m2) ^ ^ r(2 - e/2) 



x+y<i «A3-./2ni-.L(nA)i/2J r(l-e) 



Jo x+y<r'"'' A^+^{i-^y+^/^' 



where A = {x + yf - \{x + yf, ° = i - 2A- 
The integral in Eq. ( [B42|) is analyzed below. 

lo lo x^yj^''^^^^^^ [^](A-|)W 
1 , S X st 

dsdt- 



JO 



's{l - f )sl+^/2tl+^/2 (s - f - f )l+^/2 



"l /•! 1 1 



JO 



's^+H^n (1 - I - |)iW2 



(B42) 



= S + Ri + R2 + 0{e), (B43) 

where the transformations of y = st and a; = s(l — t) have been used, and the calculations for 

5", -Ri and R2 will be presented successively. 

1 ri 1 
4 



S = / dsdt- 



^0 — 
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' ds^—[- ln(- - -) + ln(l - -)] 
M 4^ ^ 4^^ 

= [/^ + 41n2(l)-/B] (B44) 
The reason for the rise of Ri and R2 can be revealed by the following decomposition: 

"1 . 1 



4 4J 

1 1 1 



/■I 1 /"^ 1 1 

= I ^^t./2(l _ |)l+e/2 X ^^iTTl[(l__±_)l+./2 - 1] 



where 



r-l 1 /.I 1 

^1 = / dt—^^ / ds 



t-/\l - |)W2 7o S^(4 - - 4Z^)^+^/2 



^^.^ Jln(4^.)-ln(3-.)] ^^^^.^4^^^^_^ (B46) 



and 



/•I 1 /■! 1 

^2= dt—T- . , , , ,„ / ds—— 

Jo t^{l- 1)1+^/2 Jo s^+^ 



= A\n{^){^)+Ic + \lA. (B47) 

The term 3 — 1.3 will be equal to 3 — 2.2, which will be shown in next two subsections. 
For the other 3-loop diagram, we can still employ the same technique. To begin with, we 
write down the corresponding integral for Fig. 15. 

g-[fc2W+(p-fe)2W](tx-tj,)g-[g^+rn2+(p-?)2+m2]txg-(p2+m2)|tj,|_ (B48) 
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Again, the time dependent part of the integral in Eq. (|B48|) can be integrated first. For 



convenience, let a = (p^ + m^), b = [{p — q)^ + m?], c = (/c^ + m^), d = [{p — kY + m^], and 
e = (p^ + m?). Then we have: 







oo 1 1 

rft^e'("+^+"+'^)*- [ (e(^+'^-'^)*- - 1) + 



c + d — e c + d + e 

(B49) 



{c + d + e){a + b + c + d)' 

In the same spirit of the calculation as performed in the previous case, we separate the 
integrals into several pieces, each of which only contains an isolated pole, and then extract the 
corresponding singular parts. 

One can start with decomposing p ■ k into: 

p.k = ^[-(p^ + + {p- kf + Sm^) + 2{p'^ + m^) + 2(P + m^)]. (B50) 
Then one can rewrite Eq. (|El8| ) as [-(3.2 - 1) + 2(3.2 - 2) + 2(3.2 - 3)], where 



3.2-1= / d'^p I d'^q [ d'^k ^'^ 



oo .y-oo J-oo {k'^ + m^){q'^ + rn?'){p^ + m?) 
1 



[fc^ + + {p — q)^ + (p — fc)^ + 4m^] 

/•oo /"OO /"OO v\ . n 

3.2-2= / d'^pl d'^q d'^k ^ ^ 



(B51) 



oo J —oo J —oo 

1 



(A;^ + rn?)[q'^ + m2)[p2 _)_ /;;2 _|_ /^^2 _|_ 3^2j 

(B52) 



+ g2 + (p*. ^2 + (p-. ^)2 + 4^2] 



-oo J-oo J-oo (^p2 _|_ ^2~)(^g2 _|_ ^2~)|p2 _|_ _|_ _|_ 3^2j 

^ (B53) 



[k^ + g2 + (p - g)2 + (p-- fc)2 + 4m2] 

Again, we have three different types of the integrals to handle 
The calculation for the term 3.2 — 1 is shown as below. 
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3.2 - 1 



d'p d'^q d'k ^ ^ 

-00 J— 00 J— 00 



1 



k"^ + + {q — py + {p — ky + Am? 

2 Xoo ^Xoo ij? + rn?){q^ + m?) Jo Xoo ^ 
1 



[/c^ — ■ k + a;p2 + xg^ — xp- q + m 



212 



2 J -00 + m2)(g2 _|_ 77^2^ 



r(2) [(a; — — ■ q + + m^]^"'^/^ 

7r'^/2r(l - e/2) .00 p.g 



-00 J ~oo 

1 



(p2 + im?){q'^ + m^) 



(X - f - |p . g + + r^)l-e/2 

7r'^/2r(l - e/2) r(2 - e/2) ^ .1 .1 1 



d'^p / d'^q / / dxdy- 
00 J-00 Jo Jo ( 



2 r(l-e/2)7-oo ^7-oo 7o ^ (g2 + m2) 

p ■ qy'"^'^ 



2 



/o i ^^^^Xl%2+L2)AW2 



{[f -(nw + i^^^r-^ 

2 io io ^2Ani-AW2(g2 + ^^2)i-. 

7r3rf/2 o ^1 ^1 l+3e/2 -e/2 

where A = x — ^ and □ = ^ — (^)^- The integration over x and y in Eq. (P54|) is carried 
out as below. 

Jo Jo '''^'^^(x|/)iW2(A-f 



'y.]l+e/2 
4 / 

X X -^1-1)1+^/2 
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Rewrite Eq. ( [B55| ) as: 

= r*(r3|^|-'"(l-J^)i; + ^ + OW]. (B56) 

One also expands the terms in Eq. ( |B55| ) up to the zero order in e. The finite part of the first 
term in Eq. ( P56|) is found to be: 



r ^?/7Y3iyP^(4 -y)- ln(3 - y)] = 4(ln + 1^- Ib- (B57) 



(1 

The second part is extracted in the same manner: 



/■' J y"^ ,f\ 1 , t, Info) ln(l-;) ,, , ^, 

A. e 



[41n(-) + -(Jc-/A)]. (B58) 



Now we turn to the term 3.2 — 2. 



3.2-2 



/oo poo roo _^ rl rl i T) • fj 

d'^p d'^q d'^k / dxdy ^ ^ 

-oo J~oc J —oo Jo Jo x+1l<:^ 



-oo J~oo J —oo JO 

1 



x+y<i 4(g2 + (g2 + 77^2) 



[A;2 — (x + y)p ■ k + {x + y)p^ + yg^ — yp- q + m^]^ 

(r(l)'^ J-00 J-oo 4(g^ + m^)Jo JO a;+s/<i r(3) 
p- q 



00 



1 f ^.^.„^r(2-e/2) r(l-e) 



00 M(g2 + m2)yo Jo x+s,<i r(3) r(2 - e/2) 
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.r(i-e) 



r(2 + e/2)r(-3e/2) 



r(l + e/2) r(l-e) 

X + J/<1 



JO "'0 x+v<i 8 



r(l + e/2) A3+^ni+^/2' 



(B59) 



Here ^ = {x + y) - \{x + yf and □ = f - 

The integration over x and ?/ in Eq. ( |B59|) is performed as below. 



X+J/<1 
1 rl 

dsdt 



A V2A' 
S X st 



-e/2^3+e 



7o (st)i+^/2A(A- f)iW2 



4 4 



4J 



(B60) 



The s dependent part in Eq. ( P60|) can be separated into: 



ds- 



1 



- 4^)1+^/2 

1 



ds 



1 
4 



(1-4^)(1 



+ / ds 



- 4^)1+^/2 

1 

4-t 



- 4^)1W2 



nl 1 

+ / ds 



s 



0(e) 



(B61) 



The first term in Eq. ( P61|) can be integrated out as: 



"'(i-i^Ki-f) 



[ln(l - 7) 



ln(l - 3)]- 



(B62) 



Inserting it into Eq. ( [B60D , we inherit 



dt(^) [ln(l -\)- ln(l - ^)] = 4(/^ -Ie) + 2 In(^) - 3 ln(^; 



(B63) 



Again, inserting the second term in Eq. (|B6 1|) into Eq. ([B6CI| ), one has: 



(B64) 
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After substituting the 3nd term in Eq. ( P61|) into Eq. ( P60|) , one obtains: 



ds- 





:— )[41n(|) 
— e 6 



dt[ 



1 - 



c 



2 

-Ia)]. 



Ml -I)] 



(B65) 



Here we evaluate the integral of 3.2 — 3. 



3.2-3 

1 



1 roo roo roo _^ 

= - d'^p d'^q / d'^k 

4 J ~oo J —CO J —CO 



p-q 



oo J-oo {p^ + m?){q'^ + rn^) Jo 
1 



dx 



{kP' — p- k + p'^ + '^"'^ 

' d^p r d^q ^- ^ 



p^ + xq^ — xp ■ q + m^)^ 

"-t J — oo J —oo 



oo J-oo {p"^ + rn?){q'^ + m?) 
1 

(fp^ — xp ■ q + a;g^ + m^)^^^/^ 

= 7 ^'^'^W^ /' r r ^2 ^^ 

4 ^0 J-oo J-oo {p^ + m?){q^ + m?) 

1 



^|)l-./2(p2 _ l^^p-. ^+ + |^2)l-./2 



"^/^ io ./o ^(g^ + m?) {p^ — ^xyp- q + ^xyq^ + ^ym?y ^1'^ 

1 ,/2r(2-6/2) r(i - 6) . . 1 



4 



(f)W2 ^'')'^'r(2-e/2)yo ^0 



dxrf?/ / d^q 



oo (g2 + m^) 



211- 



\xy)y 



y q - {\xyq) 
\{\xy - (|xy)2)g2 + 

1 (vrmi-e) , . 

4 (|)W2 y(g2 + ^2)nl-e(^2 + g^2) 

1 (7r)'^r(l - e) 2 (7r)'^/^r(2 + e/2)r(-|6 
4 (|)W2 3 r(l-e)r(l + e/2) 

' C dxdy'Ll^{^)\^. 



3n 



1-e 



(B66) 



Here □ = \xy — The integration over x and y in Eq. ( [B66|) is obtained as: 



/' ['dxdyi^r/'- 
Jo Jo 3 



xy 



^(^^)_4(^^)2]l+e/2 
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3^ Jo Jo ^^(x2/)lW2(|)l+e/2 (1 _ M)l+./2 
pi l-V 

-3/ dy[ ^ 

Jo y 

-3lE + 0{e). (B67) 



Finally, we summarize the results of two 3-loop calculations. One should observe that there 
are not any sub-divergences in 3-loop diagrams and therefore no such a term as In(cm^a^)^ 
exists. In the following summation of both 3-loop calculations results, we demonstrate this 
observation by explicit calculations. By collecting all previous results, the contributions of 
FD.311 and FD.312 in are given by: 

FD.3U:-l,ir(^„-i) Ml, -4r(^)(-i, .„(!,- 

2^r(^)41n(^)^] = (B68) 
FD.312 : -\n^HHl){-\) + 2(i)r(^)41n(^)(-l) = (B69) 

These confirm our observation mentioned above. The consequent results of the leading 
divergences which contributes to Zg are used in Eq. ^ 

APPENDIX C: BASIC INTEGRALS 

We define s) = E^=i and s) = E^=i(-l)"t^. 



Is. lo' du'^ = 2 In^d) - 4[H(|, 2) - 2)] 



Ic /o'rf«(^ = -42(1,2) 
/o^dni^ = -H(l,2) 
I,. J^du'-^ = -Eil2) 
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I^. dudv'-^^ = -[H(i, 2) + 31n(|) - 1] 
Ig. Io' du\n{u) ln(l - f ) = 2 - 31n(|) - 4S(i, 2) 

I J. /o^(i«ln(l-H) = 31n(|)-l 
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Figures Caption 



Fig. 1 

A two-dimensional cut (along a lattice plane perpendicular to the disordered substrate) of 
the three-dimensional system. 
Fig. 2 

The Feynnman diagram for the correlation function and response function. 
Fig. 3 

The Feynman diagram representing 0(x , t)(p{x, t') cos(0(x, t) — (f){x, t')). 
Fig. 4 

The basic diagram for A. 
Fig. 5 

The Feynman diagram for r2,o up to order A^. 
Fig. 6 

Two mutually canceled two-loop Feynman diagrams. 

Fig. 7 Two mutually canceled two-loop Feynman diagrams. 

Fig. 8 

A 2-loop Feynman diagram: FD. 211 
Fig. 9 

A 2-loop Feynman diagram: FD. 212 

Fig. 10 A 2-loop Feynman diagram: FD. 213 

Fig. 11 

A 2-loop Feynman diagram: FD. 214 
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Fig. 12 

A 2-loop Feynman diagram: FD. 215 
Fig. 13 

A 2-loop Feynman diagram: FD. 216 
Fig. 14 

A 3-loop Feynman diagram: FD. 311 
Fig. 15 

A 3-loop Feynman diagram: FD. 312 
Fig. 16 

A Feynman diagram: FD. nla 
Fig. 17 

A Feynman diagram: FD. nib 
Fig. 18 

Two Feynman diagrams contributing to ri^2- 
Fig. 19 

The other two Feynman diagrams contributing to Fi 2. 
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